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The most fundamental rso^irements of a radiation treatment optimization wd: (£) doss is Siomcgene- 
ously deposited in the Planning Target Volume (FTP); (ti) the dose deposited in any Organ At Bisk 
(OAB) does not exceed a threshold value and ideaffly should be zero; (ssj) the dose deposited in All The 
Kest (ATR), I.e. organs esd tissue not included in the PTV and OASte, should he as asaall as possible 
and ideally zero to minimize the risk of Gscondosy carcinogenesis; («/) the dose gradient crcaaing the 
PTIT boundaries should be es high as possible. Optmikatsons are pursmd by the wifaimiMrtioa of a 
ppsitiv<3-deGnifca objective function. A succssenxl optimization will yield a global minimum to this objec- 
tive fun&c&ora in a abort computation fcime with physically achievable beamlet intensities. 

This work presents a new approach fe> radiation treatment optimization that is very feet and yields a 
global minimiam of the objective Amotion without the use of a search routine but rather solves a linear 
system of equations. The possibility of such a direct optimization echeme has been Zmown for decades 
but it has been impossible to implement rigorously using convention®! quadratic objective functions 
because optimized results are only achievable with unphysscal negative beamlet intensities [3]. Once 
an ad-hoc condition requiring the beamlet intensities to be positive is introduced (ie. force neg&txva 
values to be zero), the linear method yields a dose distribution with artefacts. Alternatively, one may 
pursue a minimum of the objective Amotion by a direct search over c)l possible positive beam intensi- 
ties, but this £& a very time-consuming approach, Xm our work 4he negative-intensities problem is dealt 
with within tfos chjecdvz function and not as an externally imposed ad-hoc requirement. 



W<a summarize now our approach. In the following we call "organ" any organ, terget volume, region of 
tissue, identifiable anatomical entity or any defined volume within the volume espcaed to radiation. 
For simplicity we will divide the set of organs indo organs that must receive a certain doss and organs 
that should receive no dose ox a dose as small m possible- A typical objective function O satisfying the 
optimization conditions stated above is of the form: 



cilcac2£3 



where the jo A are importance ccsffidenta and tme objectivity terms ere; 



5»\A * ) 



where to,- is the weight of beamlet s, d t is the dose deposited at destination point x by beamlet i and 
is the doss r^qusred in organ h. 

The main reason to the traditional appearance of negative weights upon optimization of the obejetive 
function O is the feet that we require satisfying two conflicting on one hand we require 

O*** » 0 and on the other we require radiation to pass through the ATR (and possibly OAIte) to reach 
the PTV\ A correct requirement on O aTR (and O&Bs) ie that O****** should be minimized and Or**** 
a&ouW 6e zero only if the weighs of all thz beamlete passing through the *W<fo$e* organs are zero. 
This requirement is satiaffed if instead of that standard 0***°* above, we use new terms of the form 
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| Defeated Description of Invention » ^ 



We hav© as well added aacthsr fcann to the objects© function that replaces the unrealistic zaro-limit 
for the beesolet weights with an equal-height limit (cylindrical symmetry) x?Mch is usually the initial 
est of weights befovs optimization. This tool is of the form 



Wish the weights normalised to 

2 22 total number of beamlefcs , 
i 

O 5 ^ is positive and its minimum is zero when w t = 1 for all 0**" provides the most powerful con- 
straint to avoid negative weights. 

With the new terms introduced above, the new objecti v© function to be ussd is of the form : 



c3 oroo23 viti 


c3l ttnpsa 




A 


II 





Although not necessary (and in gsraesel count^prcdluctive) ths tesTms Qg^** can also be added to the 
objective function. In our c^caxktiosio, the terras 0^~^ were given non-zero importance cosfifcx&nts 
when we wanted to demonstrate the appearance of negative weights upon optimization. With this con- 
sideration a more general objective function is: 







fcqtsarcd <2sa# required da^? 










0000 nowel idea*** 10 



Where the coefficients j> 0 are the important parameters for organ 9 within an objectivity component 
andp^ is the importance parameter of the asymptotic symmetry condition. The optimization problem 
for all the beam intensities is reduced to the solution of a lines? system of equations as is shown in the 



The optimum of the objective function is obtained by mnis&izing the objective* Amotion Scorn the sys- 
tem of equations: 

The optimum of the positive-definite objective function O is obtained by jooinimiaing O with respect to 
all the weights w A . We perform this minimisation by requiring the set of first derivatives of O to satisfy; 

0=-^2- for all w t . 

Consider the term O^ 2 : 
dO**> 
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| Detailed Description of Invention 



In our novel approach, a&SBBftgHgB thzjonfet of th# summation to obtain; 





This simple exchange ia the summation order simplifies enormously the problem! We are aoT7 able to 
formulate the problem only in terms of the b samlet weights. The js-depsndoncs (destination point de- 

th&j*4ndep&»dent aswe: 






and 






The optimization problem for the optam 
gebraic equations: 


si weights k>, is given by the solution to ths system of linear al~ 






where 




> * 


and 






The solution to the optimisation problen 


l ia obtained by the numerical inversion of tbe matsfe a-: 




^^noud idea?*®* 1 
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§. P. GoldsasQ 0 , J. Oaea 1 aefi X jr. IBatfiSste]' 

"Dept. of Physics & Asuronomy, Uitiversity or Western Ontario, London. Ontario, Canada 

F Dept of Oncology, University of western Ontario and London Regional Cancer Centre, London, Ontsno Canada 

Afostoract 

A fast optimisation algorilhm is very important for inverse planning of Intensity Modulated Iteduuion Thenzpy (TMSCT), and for 
*»ap»ve redfoaenmy of the future. coavettUonal numerical search algorisms such as the conjugate gradient search, conducted 
with posMve beam weight constraints, generally require many iterations and may jacduce saboptimsl resufes duo so tapping in 
looaJ^iiuma A direct solution of me inverse problem using oogiventioisal quadratic objective mictions without positive beam 
C0I S ramtS ^ ?^ cfficient ^ r6SUlt m unrealistic negative beam weights. We present here a direct solution of iha inverse 
problem which does not result in tmaccepaabla negative beam weights. The objective function for the optimization of beam 
nrtmstfies for large number of beamlet* is reformulated such that the on^nizstion problem is reduced to a linear set of equations, 
ine optimal set of intensities is found through a razfrix iaversion, and negative beamlet intensities are avoided without the need for 
^^^^T^ COaS ^-?' V* m ^ cd to **** WKod to a test phantom and to a few clinical cases. We were able to achieve 
totrifeuions with very short optimization times. T^I^ optbntoion times for a single eimomical slice 
tomg a sm^e pnocesw ©>sktop computer are: 0.2 sec. for 400 beamlets; 0 sec for 1,000 beamlets; 40 sec for 2,000 beamlets and 
^min for 3,ooo beamlets. These times can be substantially fiirther isnproved using a better cjpthnizstion routine for matrct 
mvasron. In conclusion, the new method provides a fast and robust technique to find a global mmimum mas yields excellent results 
xor me inverse planning of IMRT. 



. Inverse plaaaing, ophmization, objective function. 

Intensity Modulated Radiation Therapy (hMRT) is becoming a 
new standard for radiotherapy. Given the better confbnnal dose 
distributions obtained through IMRT and its dynamic delivery 
features, adaptive radiotherapy becomes aa important fector to 
be considered. A fast optimization algorithm is crucial no* only 
for designing good radiation treatment pirns hm also for the 
successful implementation of future interactive adaptive 
treatment techniques. Conventional optimization algorithms 
using numerical searches such as the conjugate gradient search 
[1-2] with positive beam weight constraints usually require 
many iterations (i.e. long computation times) and may result in 
suboptimal plans due to snapping in local minima of the 
objective function. A direct solution of the inverse problem 
using conventional quadratic objective functions without 
imposing positive beam constraints will be computationally 
fester but will result in unrealistic negative beam weights. We 
present here a very fast method for the direct solution of the 
invesse problem (FIDO) that avoids the difficulty of negative 
beam weights and preserves efficiency. 

R5to(J8a<n>dl 

The most fundamental requirements of a radiation treatment 
ptunizatioa are: (/) dose is homogeneously deposited in the 
Planning Target Volume (PTV); (,,) the dose deposited in any 
Organ At Risk (OAR) does not exceed a threshold value and 



ideally should be zero; («D the dose deposited in All The Rest 
(ATR), Le. organs end tissue not included in the PTV and 
OARs, should ba as small as possible and ideally zero to 
minimize the risk of secondary carcinogenesis; (fv) the dose 
gradient crossing the PTV boundaries should be as high as 
possible. Optimizations ere pursued by the minimization of a 
posnivo-deftnfte objective fraction. A successful optimization 
will yield a global minimum to mis objective function in a 
short computation timo with physically achievable bsamlet 
intensities. 

This work presents a new z$pzoach to radiation seatment 
optimization that is very fast said yields a global imnimmn of 
the objective function without the use of a search routine but 
rather solves a linear system of equations. The possibility of 
such a direct optimization scheme has been known for decades 
but it has been impossible to implement rigorously using 
conventional quadratic objective functions because optimized 
results are only achievable with unphysical negative beamlet 
intensities [3]. Once an ad-hoc condition requiring the beamier 
intensities to be positive is inirodnced (i.e. force negative 
values to be zero), the linear method yields a dose distribution 
with artefacts. Alternatively, one may pursue a minimum of the 
objective function by a direct search over all possible positive 
beam intensities, but this is a very time-consuming approach. 
In our work the negatrve-intensiiies problem is dealt with 
within the objective Junction end no? as an externally imposed 
ad-hoc requirement. 
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We summarize now our approach. For simplicity we will 
consider a single PTV, a single OAR and a single ATR. A 
typical objective function O satisfying the optimization 
conditions staled above is fmerbrm; 

O « PpwOpTv + PoasPom + PatiPaxb 
where the/?* are importance ccsSicieats and the objectivity 
terms arc: 

/cQbssxalcta \' 



/cfltoaedcto V 3 



and O^- £ X "**W| . 

where u>< is the weight of beamlet i, d t is the dose deposited 
at destination point or by beamlet /end a*** is me dose 
prescribed to the PTV. The main reason for the traditional 
appearance of negative weights upon optimization of the 
obejetive function O is the feet that we inquire satisfying two 
conflicting demands: on one hand we require Oatr "0 and on 
tfie other we require radiation to pass tfcjougji mo ATR (and 
possibly OARs) to reach the PTV. A correct requirement oa 
Oath is that Oath should he minimized and Oatu should be zero 
only if the weights of all the beamtets passing through the ATR 
are zero. This requirement is satisfied if instead of 0#m wo use 
anew ATR tenn of the form 



MtATa i 

Similarly for the OAR we use: 



£ Z "WW. 



We have as well added another term to the objective function 
that replaces the unrealistic zero-limit for the baamlet weights 
with an equal-weight limit (cylindrical symmetry) which is 
usually die initial set of weights beftr* optimization. This term 
is f the form 



With the weights normalized to 

2 «&= total DumW of beamlete, 

positive and its minimum is zero whea w t »1 for all i. 
With these mediffestions, she optimization problem for all she 
beam intensities is reduced again to the solution of a linear 
system of equations of the form; 

lL a u w )~Pt (1) 

where Wj is the (unknown) weight or intensity of besmlet j, fit 
is a vector of coefficients that depends on the dose deposited 
by beam / within me PTV, and Otj is a matrix mat describes the 
product of the doses deposited by the intersecting pairs of 
beamlets i end / on different organs (each organ with its 
importance coefficient). The sez of optimal beam weights is 
obtained from (1) by inversion: 



In other words, the solution to the (Jorge) system of linear 
equations (1) is obtained quickly and accurately by inverting 
the agay cfy using standard msftjg inversion routines 



Below we present several sets of preliminary results obtained 
with this technique for a prostate case, a head and neck case, 
and an ' "interlocked rectangles" phantom on a 3D slice. Only 
primary KERMA has been included in the calculations 
reported here* but very similar results are obtained in dose- 
distributions, when the calculated team intensities sxs imported 
to a commercial treatment planning system (Therapkn Plus 
V3.S, Nucletrozi). The only difference is a a lightly diffused 
dose to the OARs due to scatter spreading effects, hi each case, 
we include calculation end optimization tunes for KBRMA 
calculations as they were obtained on a smgle-piccesscr 
desktop PC. The program was written in CM using the 
Microsoft .NET environment The msSribt inversion routine was 
obtained from Numerical Recipes (4] end translated from 
FORTRAN into C#. No effort has been devoted to maximize 
the speed of the matrix inversion procedure. 

In all cases me number of gasusy eagles used is evenly 
cSstribaSed over a Ml 360 degree circle arctand the isocestre. 
Each beam is evenly divided into beamlets of the specified 
widm resulting in the total number of beamlets quoted. The 
source to axis distance (SAD) is 100 cm. m each case we also 
present fee DVH, scaled to 100% volume oa the vertical axis 
and 100% dose on the horiaoa&sfl asus, and a colour-coded dose 
deposition map with blue shewing no dose deposition (primary 
KERMA only) end red the largest dose deposition. 



We have developed a last and robust technique co find a global 
minimum mat yields excellent results for the inverse 
optimization problem for the rsdiatiea treatment of turnout, , 
using 1 



Work is currently in progress on a fid! inmlemeirt^on HDO in 
our treatment planning system that uses collapsed cone 
convolution method for dose calculation [5]. Work is 
proceeding as well on an efficient 3D implementation for on- 
line adaptive radiotherapy as might be possible with helical 
rancmerapy [6]. 
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Research and Development Challenge Fund (OCITS project). 
We also would like to thank Dr. E. Wong for useful 
discussions and for suggesting the challenging geometsy of the 
interlocked rectangle phantom. 
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Figure 1: Head and Neck Cose. - 20 gantry angles 




Figure 3i Interlocked Rectangle Phantom. - 20 gantry 
angles 




Figure 2: Head and Neck Case. - 40 gantry angles 
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Set No. 4: Interlocked Rectangles. -40 gantry angles 
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Set No. 5: Prostate Case (Panoramic and dose-up views). 
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Wo claim: 

1. A method for optimizing radiation treatment, said method comprising the 
step of resolving an objective function, said objective function comprising 
importance coefficients and objectivity terms. 

2. A system for optimizing radiation treatment, said system comprising 
means for resolving an objective function. 

3. A method of planning delivery of radiation therapy to maximize radiation to 
a pfanned target volume and minimize radiation to surrounding tissues 
outside the planned target volume, said method comprising the step of 
resolving an objective function. 

4. The invention substantially as described and illustrated herein. 
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